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SECTION  1 
BACKGROUND 


In  the  source  region  of  a  nuclear  burst  an  electric  field  is 

generated  by  the  outward  flux  of  Compton  electrons.  This  field  is  limited 

by  the  conductivity  of  the  air  which  is  a  result  of  the  ionization  produced 
by  the  Compton  electrons.  The  field  builds  up  until  the  conduction 

electrons  drift  inward  at  a  rate  that  compensates  the  outward  displaced 

Compton  electrons.  This  equilibrium  field  is  given  by 

E  =  J  /a  (1) 

s  c 

where  Jc  is  the  Compton  current  and  a  is  the  conductivity  of  the  air. 

The  conductivity,  when  the  motion  of  the  ions  are  neglected,  is 

given  by 

a  =  n  e  u  (2) 

where  n  is  the  free  electron  density  and  the  electron  mobility,  u,  is  a 
function  of  E.  In  general,  the  electron  density  depends  on  the  dose  rate, 
y,  the  electron  attachment  rate,  a,  the  ion-ion  recombination  coefficient, 
0.,  as  well  as  the  electron-ion  recombination  coefficient,  e  .  These  para¬ 
meters  are  functions  of  electric  field  and  the  water  vapor  content  of  the 


At  low  y,  n  is  determined  by  the  attachment  rate.  In  steady  state 
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where  K  =  2.2  x  10 9  cm-3/rad  is  the  generation  constant  in  air.  Since  the 
dose  rate  is  proportional  to  the  Compton  current. 


Jc  =  Kc  * 


The  equilibrium  field  is  then 


which  is  independent  of  y  and  is  called  the  saturated  field. 

At  high  y,  n  is  determined  by  electron-ion  recombination.  In 
steady  state 


■  ■ 


The  equilibrium  field  is  now 


s  e 


Since  this  field  increases  with  y1^2,  it  is  not  truly  saturated.  The  high 
dose  rate  regime  is  defined  by 


n  e  »  a 
e 


An  initial  value  of  “2  x  10-7  cm3/s,  which  is  due  to  two-body 
recombination  with  N2+  and  02+,  has  been  measured  at  low  pressure  mainly 
using  drift  tubes.  The  value  of  increases  during  the  pulse  as  ions  form 
agglomerates  like  N3+,  N4+  as  well  as  clusters  with  water  molecules.  In  a 
typical  calculation  will  reach  a  plateau  value  of  »1  xl0~6  cm3/s  for 
moist  air  (Ref.  5).  Recently  there  has  been  theoretical  and  experimental 
work  on  another  mechanism  which  is  possible  when  complex  molecules  like  H20, 
C02,  NH3,  etc.,  are  present  in  large  densities  to  act  as  a  third  body  in  the 
recombination  process. 

1.1  TWO-  AND  THREE-BODY  RECOMBINATION. 

The  free  electrons  in  a  gas  follow  random  paths  caused  by  elastic 
collisions  with  the  gas  molecules.  If  an  electron  approaches  to  within  some 
sphere  around  a  positive  ion  there  will  be  an  attractive  force  between  them. 
The  radius  of  this  sphere  is  usually  taken  to  be  when  the  electrostatic 
energy  is  equal  to  the  average  kinetic  energy  of  the  electrons,  that  is 


Recombination  may  still  not  occur.  The  probability  of  various  recombination 
mechanisms  (radiative,  dissociative,  or  dielectronic)  are  given  by  the 
quantum  mechanical  cross  sections.  The  two-body  process  can  be  represented 


where  B2  1S  the  two-body  rate  constant.  Increasing  the  gas  density  does  not 
affect  g2  since  the  motion  is  already  random,  and  a  collision  is  as  likely 
to  lead  to  the  escape  of  the  electron  as  to  the  capture. 

If  electrons  can  undergo  inelastic  collisions  with  neutral  gas 
molecules,  then  three-body  recombination  can  taken  place.  This  can  be  rep¬ 
resented  as 


&3 

e  +  A+  +  M  - -  A  +  M*  (11) 

where  B3  is  the  three-body  rate  constant.  At  low  pressure  the  rate  of  this 
process  is  directly  proportional  to  the  density  of  the  third  body  because 
the  probability  of  recombination  is  greatly  increased  if  the  electron  loses 
a  large  fraction  of  its  kinetic  energy  while  at  a  radius  less  than  r0.  At 
high  pressure  (often  near  1  atm)  the  process  saturates  and  recombination 
becomes  controlled  by  the  electron  mobility  (or  diffusion  constant).  The 
total  recombination  coefficient  (below  saturation  by  the  third  body)  is 

3e  =  e2  +  83  N  (12) 

wnere  N  is  the  number  density  of  the  third  body.  When  the  third  body  is  a 
non-polar  molecule  like  C02,  vibrational  excitation  is  the  important 
inelastic  process,  but  for  polar  molecules  like  H20,  rotational  excitation 
dominates. 


1.2  PREVIOUS  MEASUREHENTS. 

Microwave  absorption  methods  have  been  used  to  measure  a  number  of 
two-  and  three-body  recombination  rates  (Ref.  8).  Measurements  in  pure 
water  vapor  up  to  20  Torr  at  21 °C  and  81 °C  (Ref.  6)  were  used  to  deduce 
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From  the  three-body  theory  it  is  not  expected  that  the  type  of  ion  would  be 
important  in  determining  e3,  so  the  value  measured  for  pure  water  should  be 
applicable  to  moist  air  using  the  H2O  number  density  for  N. 

1.3  FIELD  DEPENDENCE. 

However,  the  measurement  of  &3  was  performed  at  essentially  zero 
electric  field,  while  the  b3  of  interest  to  EMP  predictions  might  be  «1  kV/ 
cm.  Fields  of  this  size  heat  the  electrons  to  a  few  tenths  of  an  eV.  The 
dependence  of  b3  on  gas  temperature  given  above  is  primarily  a  reflection  of 
the  dependence  of  e3  on  the  electron  temperature.  The  three-body  rate  con¬ 
stant  is  expected  to  go  as  rQ3  (Ref.  6),  and  Equation  (9)  shows  that  rQ3 
decreases  as  the  temperature  of  the  electrons  increases.  Let  us  assume  that 
the  T-2  dependence,  derived  by  changing  T  by  21%,  holds  over  an  electron 
temperature  increase  of  a  factor  of  four  to  0.1  eV.  Then  b3  would  be  at 
least  a  factor  of  16  smaller  than  that  found  above  and  Be  would  not  be 
much  larger  than  the  two  body  rate. 

The  only  previous  measurement  of  a  three-body  electron-ion  recom¬ 
bination  rate  made  at  a  non-zero  field  was  performed  on  C02  (Ref.  4).  Here 
it  was  found  that  the  three-body  rate  fell  very  rapidly  with  electric  field 
while  the  two-body  rate  fell  only  very  slowly. 


SECTION  2 
EXPERIMENT  DESIGN 


Although  the  dose  to  the  air  in  a  nuclear  blast  is  initiated  by 
the  prompt  and  delayed  gamma  rays,  it  is  the  Compton  electrons  that  produce 
the  free  electron-ion  pairs  that  cause  the  air  to  be  conducting.  In  order 
to  reach  the  high  dose  rate  regime,  where  recombination  is  the  dominant  free 
electron  removal  mechanism,  in  an  above  ground  test,  it  is  necessary  to  use 
an  electron  beam  directly  from  a  simulator. 

As  we  have  seen,  the  measurement  of  the  electron-ion  recombination 
coefficient  must  be  performed  at  realistic  fields  in  order  to  be  relevant  to 
EMP.  This  can  be  done  with  a  parallel -pi ate  ion  chamber.  The  essential 
idea  of  the  experiment  then  becomes: 

1.  Apply  an  electric  field  to  the  gas  in  the  ion  chamber. 

2.  Shoot  a  pulse  of  high  energy  electrons  through  the  chamber  to 
produce  a  large  density  of  electron-ion  pairs. 

3.  Measure  the  ion  chamber  current  during  the  pulse  to  infer  the 
electron  density. 

4.  Calculate  ee  using  previously  measured  air  properties 
(e.g.,  electron  mobility  and  attachment). 
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2.1  THE  ELECTRON  BEAM. 


HIFX  (flash  x-ray  machine  at  Harry  Diamond  Laboratories)  was  used 
as  the  electron  beam  generator.  With  a  4.1  MV  charge,  HIFX  produces  a  fairly 
narrow  distribution  of  electron  energies  with  a  peak  at  2.4  MeV  and  a  maxi¬ 
mum  energy  of  3.0  MeV.  The  beam  is  transported  in  a  6"  diameter  pipe  of 
variable  length  and  pressure.  Extremely  high  dose  rates  can  occur  when  con¬ 
ditions  are  such  as  to  produce  a  beam  pinch.  Under  more  stable  conditions 
it  is  possible  to  produce  dose  rates  of  10 14  rad/s.  Most  of  the  data 
reported  here  is  for  a  dose  rate  of  4  x  10 13  to  10 14  rad/s  at  the  entrance 
to  the  ion  chamber.  At  HIFX  energies,  the  value  of  Kc  in  Equation  (4)  is 
5.9  x  10- 12  (A/cm2) /( rad/s) ,  so  beam  current  density  was  as  large  as 
600  A/ cm2. 

If  we  use  a  nominal  value  of  a  =  108  s_1  and  a  minimum  value  of 
2  x  1 0-6  cm3/s  for  0e  in  Equation  (8)  we  find  that  recombination  dominates 
over  attachment  for  dose  rates  large  compared  to  2  x  10 12  rad/s.  Thus  the 
HIFX  dose  rates  are  quite  adequate  for  studying  recombination. 

The  shape  of  the  HIFX  beam  current  waveform  varied  extensively 
from  shot  to  shot.  The  FWHM  varied  from  10  to  20  ns  but  was  usually  around 
16  ns.  The  rise  time  was  also  variable  and  the  pulse  often  rose  in  two  or 
more  steps.  At  a  dose  rate  of  5  x  10 13  rad/s  for  a  time  of  10  ns,  10 15  cm- 3 
free  electrons  are  generated  in  the  air.  Later  it  will  be  shown  that  the 
equilibrium  electron  densities  were  on  the  order  of  1014  cm-3,  so  that  by 
the  peak  of  the  HIFX  pulse  there  is  sufficient  time  for  the  electron  density 
to  grow  to  the  point  that  it  is  controlled  by  recombination. 


2.1.1 


Beam  Transport. 


If  a  high  current  electron  beam  is  injected  into  a  vacuum,  the 
strong  electrostatic  repulsion  will  cause  the  beam  to  "blow  up".  Only  if  a 
gas  is  used  to  produce  a  background  of  positive  ions  can  the  repulsive  force 
be  reduced  enough  that  it  is  balanced  by  the  inward  force  of  the  beam's 
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self-magnetic  field.  Such  a  "charge  neutralized"  beam  can  then  be  propa¬ 
gated  down  a  conducting  drift  tube  with  the  return  current  flowing  on  the 
tube's  inner  surface.  If  the  gas  density  is  too  high,  the  magnetic  force 
will  cause  the  beam  to  constrict  and  a  pinch  may  occur. 

In  some  electron  beams  (especially  those  of  large  v/y,  see  below) 
there  is  also  a  tendency  for  the  beam  to  be  "current  neutralized".  Because 
of  the  changing  flux  associated  with  the  beam's  current  the  return  current 
tends  to  flow  as  close  to  the  beam  current  as  possible.  When  the  conduc¬ 
tivity  of  the  gas  is  large  enough  most  of  the  return  current  can  be  carried 
by  the  gas,  reducing  the  inductance  of  the  beam  plus  return  path. 

A  critical  parameter  in  the  transport  of  intense,  relativistic 
electron  beams  is  called  v/y  (Ref.  2).  y  is  the  usual  relativistic  measure 
of  a  particle's  energy  in  units  of  its  rest  energy,  v  measures  the  "linear 
density"  of  a  beam  and,  as  such,  is  a  measure  of  the  electrostatic  repul¬ 
sion.  v  is  defined  as  the  number  of  electrons  per  cm3  times  the  classical 
electron  radius.  The  ratio  v/y  is  then  a  measure  of  the  "stiffness"  of  the 
beam;  the  smaller  v/y,  the  stiffer  the  beam  and  the  easier  to  get  the  beam 
to  propagate.  It  has  been  shown  that  v/y  must  be  less  than  one  for  a  beam 
to  propagate  without  an  external  magnetic  field. 


The  HIFX  beam  has  y  =  5.7.  The  maximum  current  output  is  20  kA. 
This  gives  a  maximum  value  of  v/y  of  0.2,  so  the  HIFX  beam  is  quite  stiff 
and  easy  to  propagate  with  the  proper  background  pressure.  When  the  beam  is 
passed  through  a  collimator,  v  is  reduced,  and  the  beam  passing  into  the  ion 
chamber  has  a  v/y  of  only  0.02. 


Two  configuations  were  used  in  this  experiment.  In  the  first,  the 
outer  chamber  that  housed  the  ion  chamber  was  mounted  to  the  end  of  the 
drift  tube,  and  the  beam  passed  directly  from  the  partial  pressure  of  the 
drift  tube,  through  a  2-mil  A1  foil  window,  into  the  collimator  and  ion 


chamber  containing  1  atm  of  air.  In  the  second  configuration,  the  chamber 
was  separated  from  the  drift  tube,  which  consisted  only  of  the  10-cm  pumping 
section.  The  1  1/4"  air  gap  between  the  2-mil  Ti  window  on  the  drift  tube 
and  the  chamber  window  was  surrounded  by  a  1  1/8"  Cu  pipe.  Four  Cu  strips 
connected  from  the  chamber  to  the  drift  tube.  Dosimetry  results  presented 
later  show  that  the  additional  air  gap  introduced  a  significant  divergence 
to  the  beam  as  it  passed  through  the  ion  chamber.  Presumably  this  was  due 
to  a  pinching  (and  then  a  spreading)  of  the  beam  in  the  air  gap  or  to  a 
perturbation  of  the  beam  return  current,  since  the  effect  was  too  large  to 
be  explained  on  the  basis  of  additional  scattering  in  the  Ti  window  and  air. 

In  the  course  of  the  experiment  it  was  found  that  the  best  pro¬ 
cedure  was  to  operate  the  drift  tube  close  to  high  vacuum  conditions  and 
vary  the  length  of  the  tube  to  vary  the  dose.  The  dose  at  the  end  of  a 
fixed  length  of  drift  tube  could  be  increased  by  increasing  the  pressure  in 
the  tube,  but  this  caused  very  sharply  peaked  pulses  which  make  the  data 
very  difficult  to  analyze. 

2.1.2  Charge  Neutral ization. 

The  charge  injected  into  the  chamber  by  the  beam  would  produce 
very  large  fields  if  the  beam  were  not  neutralized.  For  example,  if  the 
beam  was  1  cm  in  radius  and  carried  500  A/cm2,  the  longitudinal  electric 
field  at  the  edge  of  a  1  cm  slice  of  the  beam  would  be 


2  ve0 

=  95  kV/cm 


This  is  also  the  value  of  the  radial  E  field  at  the  edge  of  the  beam. 

If  the  electron  density  in  the  gas  is  10 14  cm-3  the  conductivity 
is  about  0.02  S/cm.  The  charge  relaxation  time  is  then 


t  *  e0/a 
=  5  ps 


(15) 


So  if  the  rise  time  of  the  pulse  is  5  ns,  the  fields  reach  only  values  close 
to  100  V/cm.  If  this  slice  of  beam  lies  between  the  center  and  ground  elec¬ 
trodes  of  the  ion  chamber,  neither  the  radial  or  longitudinal  fields  are 
directed  such  as  to  drive  current  across  the  gap  of  the  ion  chamber.  The 
effect  of  the  charge  carried  by  the  beam  was  observed  in  the  voltage  measur¬ 
ing  grids  in  the  chamber  and  will  be  discussed  further  under  the  data  analy¬ 
sis. 


2.1.3  Back  emf  of  the  Beam. 

As  mentioned  above,  the  circuit  consisting  of  the  beam  and  return 
current  path  has  an  inductance.  A  change  in  the  beam  current  produces  a 
back,  longitudinal  electric  field  that  does  drive  current  in  the  air  between 
the  plates  of  the  ion  chamber.  The  inductance  has  two  components:  one  due 
to  magnetic  flux  in  the  beam  and  one  due  to  the  flux  that  lies  between  the 
beam  and  the  return  path.  The  current  density  is  approximately  constant 
over  the  cross  section  of  the  beam.  The  back  emf  varies  across  the  beam, 
being  a  maximum  at  the  center  and  zero  at  the  surface  of  the  return  path. 

In  Figure  1  the  dashed  line  shows  the  path  to  which  we  apply 
Faraday's  Law. 
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Figure  1.  A  back  emf  is  produced  by  a  change  in  the  current 
of  the  e-beam. 

This  is  the  maximum  field  which  occurs  on  the  axis  of  the  beam.  When  aver¬ 
aged  over  the  cross-sectional  area  of  the  beam,  the  field  is 


E  = 


(17) 


In  order  to  keep  E  small,  a  should  be  kept  small.  That  is,  the  chamber  must 
be  designed  so  that  the  beam  return  current  flows  as  close  to  the  beam  as 
possible.  In  our  case,  r0  =  1  cm,  a  =  3  cm,  I  »  10 11  A/s.  This  gives  E  = 
140  V/cm.  This  is  a  small  effect  at  fields  of  prime  interest  to  EMP 
(1000  V/cm)  but  could  affect  measurement  at  lower  applied  fields.  The 
effect  of  this  field  is  reduced  because  the  ion  chamber  used  here  is  a 
double-sided  chamber.  Since  the  back  field  is  direct  along  the  beam,  in  one 
half  of  the  chamber  it  adds  to  the  applied  field  while  in  the  other  half  it 
subtracts.  A  smaller  compensating  effect  is  due  to  the  fact  that  as  the 
electric  field  increases,  the  electron  mobility  decreases. 


The  back  emf  that  has  just  been  calculated  is  the  force  that 
causes  electron  beams  to  become  current  neutralized.  The  backward  driven 
plasma  current  in  the  ion  chamber  is 
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J  =  a  E 
=  3  A/ an 2 

Since  this  is  negligible  compared  to  the  beam  current  density,  the  HIFX  beam 
is  not  current  neutralized  in  the  ion  chamber,  and  a  Rogowski  coil  can  be 
used  to  measure  the  beam  current. 

2.2  THE  CHAMBER. 

An  outer  chamber  was  made  to  serve  as  the  vacuum  container  for 
either  of  two  ion  chambers  with  different  electrode  spacings  or  an  equilib¬ 
rium  field  diode  (to  be  discussed  below).  A  scale  drawing  of  this  chamber 
is  shown  in  Figure  2. 

The  chamber  was  made  to  seal  to  the  end  of  the  drift  tube.  A  1/2" 
A1  plate  backed  by  a  1/8"  Pb  sheet  was  mounted  on  the  front  of  the  chamber 
with  A1  bolts  to  serve  as  a  stopping  block  for  the  beam  that  did  not  enter 
the  chamber.  An  internal  collimator,  with  a  3/4"  diameter  hole  consisted  of 
two  pieces:  one  also  acted  as  a  clamp  ring  to  hold  the  front  window  and  the 
other  as  a  housing  for  the  Rogowski  coil.  The  thickness  of  the  Rogowski 
coil  housing  was  varied  to  accomodate  the  different  thickness  ion  chambers. 

The  chamber  was  vacuum  tight  and  bakeable  to  greater  than  200°C  in 
order  to  remove  all  water  vapor  before  filling  with  dry  air.  Due  considera¬ 
tion  was  taken  of  thermal  expansion.  In  particular,  the  outer  area  of  the 
backplate  was  made  1/8"  thick  in  order  to  flex  at  high  temperature. 

There  were  two  2-mil  A1  foil  windows.  The  front  window  separated 
the  partial  pressure  of  the  drift  tube  from  the  vacuum  or  atmospheric  pres¬ 
sure  of  the  chamber.  The  window  behind  the  ion  chamber  separated  the  cham¬ 
ber  air  from  the  beam  diagnostics  which  were  located  in  the  room  air.  Both 
windows  were  held  in  place  by  clamping  the  foils  against  Viton  0-rings. 
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There  was  one  high  voltage  feedthrough  and  six  SMA  vacuum  feedthroughs 
mounted  in  the  backplate. 


The  signals  were  carried  on  Cujac  cables  from  a  junction  box 
mounted  on  the  rear  of  the  chamber  to  a  second  junction  box  where  they  were 
joined  to  RG-223  cables.  These  cables  ran  in  high  quality  (low  transfer 
impedance)  conduit  to  the  screen  room.  The  Cujac  cables  were  shielded  from 
stray  bremsstrahlung  radiation  with  a  layer  of  lead  bricks.  The  conduit  was 
connected  to  the  screen  room  and  then  carried  inside  and  connected  to  a 

patch  box  in  an  equipment  rack.  The  cables  in  the  conduit  connected  to 

panel  feedthroughs  on  the  front  of  the  box,  and  RG-58  cables  were  then  con¬ 
nected  from  the  patch  panel  to  the  various  recording  devices.  The  most 
important  of  the  data  were  recorded  on  a  Tektronix  7912  AD  transient  digiti¬ 
zer,  often  with  two  signals  recorded  sequentially.  The  other  fast  signals 
were  recorded  on  Tektronix  7844  oscilloscopes. 

2.2.1  Air  Samples. 


The  chamber  was  baked  at  120°C  for  1  hour  while  under  diffusion 
pump  vacuum.  When  the  temperature  returned  to  that  of  the  room,  the  pres¬ 
sure  read  2  x  10- 5  Torr.  The  chamber  was  then  filled  with  commercial  dry 
air  to  a  pressure  of  1  atm.  Although  the  chamber  was  originally  designed  to 


measure  the  water  content  of  the  air  from  IR  absorption,  the  wet  air  data 
was  obtained  on  air  saturated  with  water  vapor.  The  chamber  was  filled  with 
water  vapor  from  a  reservoir  and  then  filled  to  1  atm  with  room  air  which 
already  had  an  80%  or  greater  relative  humidity.  The  excess  water  condensed 
on  the  internal  chamber  surfaces,  and  the  air  remained  saturated  at  the 
chamber's  temperature  during  the  experiment.  The  chamber  temperature  was 
measured  with  a  calibrated  thermistor. 


:vw 


2.2.2 


Ion  Chambers. 


The  ion  chamber  is  shown  schematically  in  Figure  3.  The  ion  cham¬ 
ber  is  double-sided  with  a  center  electrode  and  two  ground  electrodes  of 
1  mil  A1  foil.  There  are  two  voltage  measuring  grids  made  from  90%  trans¬ 
parent,  5  mil,  A1  wire  mesh.  The  electrodes  and  grids  were  separated  with 
teflon  spacers.  Two  spacings  were  used:  3/16"  for  the  "1/2  cm"  chamber  and 
3/8"  for  the  "1  cm"  chamber.  In  order  to  reduce  the  inductance,  the  chamber 
measuring  resistor  was  made  up  of  about  40  parallel  resistor  placed  radially 
outward  from  the  center  electrode.  Also  the  storage  capacitor  was  a  large 
number  of  parallel  capacitors  placed  around  the  periphery  of  the  ion  cham¬ 
ber.  Both  chambers  had  a  capacitance  of  0.032  yF.  The  "1/2  cm"  chamber  had 
a  resistance  of  3.40  n  and  the  "1  cm"  chamber  had  a  resistance  of  3.18  a. 

The  ion  chamber  capacitors  also  serve  as  the  return  path  for  the 
beam  current.  During  the  course  of  the  experiment  an  oscillation  occurred 
in  the  ion  chamber  current  for  very  fast  rising  pulses.  Ten-ohm  resistors 
were  inserted  into  the  ground  leads  of  one  set  of  capacitors,  in  order  to 
put  some  damping  resistance  into  the  circuit  consisting  of  beam  inductance 
and  storage  capacitance.  Also  a  Cu  foil  was  wrapped  over  the  outside  of  the 
capacitors  to  confine  any  leakage  flux.  This  served  to  effectively  remove 
the  oscillations. 

2.2.3  Scattering. 

The  electrons  passing  through  the  foils  and  grids  of  the  ion  cham¬ 
ber  undergo  multiple  Coulomb  scattering.  The  divergence  of  the  beam  is 
increased  as  it  passes  through  each  window  and  electrode.  If  we  assume  that 
there  is  effectively  4  mil  of  A1  at  the  front  of  the  chamber,  then  the  rms 
scattering  angle  is  8°,  and  the  beam  area  is  increased  by  about  30%  at  the 
rear  of  the  chamber.  The  divergence  of  the  beam  was  measured  and  accounted 
for  in  the  data  analysis. 
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2.2.4  Boundary  Layer. 


Very  large  densities  of  electrons  and  ions  are  created  in  this 
experiment.  When  the  electrons  drift  under  the  motion  of  the  applied  field, 
a  layer  of  positive  ions  are  left  behind  at  the  cathode.  The  potential  can 
quickly  appear  across  the  boundary  layer  rather  than  the  bulk  of  the  gas. 
As  an  example  let  us  assume  that  the  ions  are  immobile,  n  =  1014  cm-3,  y  = 
103  cm2/V*s,  and  that  103  V  are  applied  to  a  1  cm  gap.  Then  in  5  ns  the 
electrons  would  leave  behind  a  layer  of  positive  ions  5  x  10- 3  cm  thick  with 
a  voltage  of  more  than  2  x  10 3  V  across  it.  Since  this  is  more  than  the 
available  voltage,  there  would  be  no  field  on  the  bulk  of  the  gas. 

In  the  above  example,  as  the  electric  field  builds  in  the  boundary 
layer,  the  ion  current  to  the  cathode  increases.  Since  the  ion  mobility  is 
about  103  times  smaller  than  the  electron  mobility,  the  field  in  the  boun¬ 
dary  layer  will  greatly  exceed  that  in  the  bulk  of  the  gas  before  the  ion 
current  to  the  cathode  can  match  the  electron  current  in  the  bulk.  Often 
the  electron  avalanche  field  will  be  exceeded  in  the  boundary  layer.  There 
are  then  a  number  of  processes  that  determine  the  boundary  layer  voltage: 

1.  Ions  in  the  boundary  layer  drift  to  the  cathode. 

2.  A  small  fraction  (=0.001)  of  the  ions  cause  electrons  to  be 
emitted  upon  impact. 


3.  The  emitted  electrons,  as  well  as  free  electrons  produced  by 
the  radiation,  avalanche  in  the  boundary  layer. 


4.  This  increases  the  ion  density,  decreases  the  thickness  of 
the  boundary  layer,  and  reduces  the  boundary  layer  voltage. 


A  detailed,  computer  calcuation  is  necessary  to  predict  the  boundary  layer 
vol tage. 
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In  this  experiment  grids  were  used  to  measure  the  electric  field 
in  the  bulk  of  the  gas.  The  center  electrode  was  usually  operated  at  a  neg¬ 
ative  voltage  and  the  voltage  from  grid  to  ground  (anode)  does  not  include 
the  voltage  on  the  boundary  layer. 

2.2.5  Grids. 


The  Debye  shielding  length  is  given  by 


(18) 


In  the  fields  applied  during  this  experiment  the  electrons  are  heated  to 
tenths  of  an  eV.  Equation  (18)  then  gives  a  Debye  length  on  the  order  of 
10-5  to  10-4  cm.  Since  this  is  much  smaller  than  any  of  the  physical  dimen¬ 
sions  of  the  ion  chamber,  the  ionized  air  in  this  experiment  must  be  treated 
as  a  plasma. 


The  voltage  measuring  grids  are  really  probes  immersed  in  the 
plasma.  The  Debye  length  is  much  smaller  than  the  spacing  between  the  grid 
wires,  and  each  wire  acts  independently.  The  constraints  on  the  grids  are 
then: 


1.  The  grids  must  be  placed  along  an  equi potential . 

2.  The  grids  must  be  sufficiently  transparent  to  not  impede  the 
current. 


3.  The  grids  must  remain  close  to  the  plasma  potential. 


i 

i 


The  grid  circuit  used  for  the  data  taken  during  this  experiment  is 
shown  in  Figure  4. 
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Figure  4.  Grid  circuit. 

Ra  is  the  resistance  of  the  air,  is  the  capacitance  of  the  grid  to 
ground,  and  -V  is  the  voltage  on  the  bulk  of  the  air.  Grid  current  is  drawn 
through  both  the  resistive  divider  and  the  stray  capacitance.  If  the  plasma 
potential  at  the  grid  is  500  V,  the  resistor  carries  100  mA.  The  stray 
capacitance  of  the  grid  to  ground  is  about  10  pF.  If  the  grid  is  required 
to  charge  to  the  500  V  in  1  ns,  then  the  current  drawn  by  the  capcitance  is 
5  A,  much  larger  than  the  resistor's  current. 


This  current  must  be  furnished  by  the  plasma.  The  random  electron 
current  available  in  the  plasma  is 


J 


e 


n  e  v 

e 


(19) 


270  A/cm2 


The  random  positive  ion  current  is  less,  both  because  the  ion  mass  is  much 
larger  and  because  the  ion  temperature  is  thermal. 


J.  =  N+  e 


(20) 


=  0.2  A/ cm 2 

The  grid  has  a  surface  area  in  the  plasma  of  2.85  cm2  x  0.1  x  rr, 
or  0.9  cm2.  The  electron  current  is  initially  much  larger  than  the  ion  cur¬ 
rent  and  much  larger  than  the  current  carried  to  ground  by  the  grid.  After 
the  grid's  capacitance  has  been  charged,  the  potential  of  the  grid  adjusts 
negatively  with  respect  to  the  plasma  to  a  value  which  cuts  the  electron 
current  such  that 

I  -  I.  =  100  mA 
e  i 

If  we  assume  that  the  electron  distribution  is  still  Maxwellian,  despite  the 
heating  by  the  field. 


-  Si 

kT 

0.3  A  =  270  A/cm2  x  0.9  cm2  x  e 

and  V  =  7  V.  So  the  grid  is  driven  only  slightly  below  the  plasma  potential 
at  the  grid's  location. 
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2.3 


THE  MEASURING  CIRCUIT. 


There  were  very  stringent  design  requirements  placed  on  the  exper¬ 
iment  by  the  measuring  circuit.  The  logical  flow  is  as  follows: 

1.  The  large  amount  of  ionization  produces  a  very  low  air  resis¬ 
tance  in  the  chamber,  =2  q. 

2.  The  measuring  resistor,  between  the  center  electrode  and  the 
storage  capacitor  must  be  on  he  order  of  the  air  resistance 
or  most  of  the  applied  voltage  appears  on  the  measuring 
resistor  instead  of  the  chamber. 

3.  But  the  measuring  resistor  also  serves  as  the  damping  resis¬ 
tor  in  the  circuit  consisting  of  the  chamber  capacitance,  the 
storage  capacitor,  and  the  wiring  inductance. 

4.  Thus  to  prevent  oscillations  the  inductance  in  the  measuring 
circuit  must  be  held  to  a  minimum. 

As  an  example,  we  estimate  the  capacitance  of  the  chamber  to  be 
10  pF  and  the  inductance  of  a  wire  running  to  a  resistor  and  storage  capaci¬ 
tor  eternal  to  the  chamber  to  be  100  nH.  The  frequency  of  oscillation  is 
then  150  MHz,  and  the  critical  damping  resistance  is  200  a.  Either  the  cir¬ 
cuit  will  oscillate  with  a  small  measuring  resistor  or  almost  no  voltage 
will  be  across  the  chamber  with  a  large  measuring  resistor.  It  is  therefore 
necessary  to  locate  the  resistor  and  capacitor  inside  the  vacuum  chamber  and 
to  distribute  them,  as  described  above,  in  order  to  reduce  the  inductance. 

If  the  inductance  was  reduced  to  0.1  nH,  then  the  critical  damping 
resistance  would  be  6  q.  A  compromise  value  of  =3  n  was  used.  The  circuit 
is  sufficiently  damped,  and  only  about  one  half  of  the  applied  voltages  is 
lost. 


2.3.1 


Circuit  Analysis. 


The  circuit  used  to  measure  the  ion  chamber  current  is  shown  in 


Ct  =  storage  capacitor  =  0.032  pF 
C2  -  coupling  capacitor  =  0.001  pF 
Cs  =  air  sample  capacitance  •  10  pF 
R,  =  measuring  resistor  =  3.4  a 
R2  =  terminating  resistor  =  50  a 
Rs  =  air  sample  resistance  •  2  a 
»s  =  sample  (chamber)  voltage 
=  measured  voltage 

Figure  5.  Measuring  circuit. 

Not  shown  on  the  diagram  is  the  voltage  source  used  to  initially  charge  all 

capacitors  to  V^;  the  power  supply  is  resistively  isolated  from  the  circuit 

on  the  time  scale  of  interest.  Also  not  shown  is  a  voltage  monitor  circuit 

used  to  verify  that  a  voltage  of  V  appeared  at  the  chamber.  V  (t)  was 

o  m 

recorded  and  we  wish  to  infer  I ( t)  =  ii  +  i2  +  i<j,  the  total  current  flowing 
through  the  air  in  the  chamber. 

Before  presenting  the  exact  equations,  it  is  of  value  to  think  of 
the  circuit  in  terms  of  three  time  constants.  A  change  in  the  air  resis¬ 
tance  causes  the  voltage  across  the  chamber  to  change  with 


Figure  5. 


where  the  ||  symbol  is  used  to  represent  the  reciprocal  of  the  sum  of  the 
reciprocals  of  two  numbers.  This  time  constant  shows  that  V ( t)  will  closely 
follow  R<Jt).  The  measured  voltage  is  related  to  V<.  by  a  time-constant  dif¬ 
ferentiation  with  a  time  constant  of 

t  2  =  R2C2  (22) 

=  50  ns 

Since  the  HIFX  pulse  had  a  FWHM  of  20  ns,  a  small,  but  not  negligible,  decay 

occurs  on  V  when  V_  is  unchanged.  The  third  time  constant  controls  the  sag 

m  i> 

of  the  applied  voltage 

t3  -  (Rx  +  R$)  x  Cx  (23) 

=  170  ns 

This  is  quite  Tong  compared  to  the  20  ns  HIFX  pulse,  so  the  voltage  applied 
to  the  resistor  combination  of  Rj  and  R<.  stays  constant  at  Vq. 

If  we  do  not  include  the  time  constant  effects,  then  the  equations 
are 


VS  = 
1 1  = 

*2  = 


V  +  V 
o  m 

-VRi 
m  1 

_V_/R2 

m 


dV 


ic  =  -C< 


m 


dt 


(24) 


The  current  i$  is  always  small  (less  than  1%)  compared  to  i  lt  especially 
since  measurements  are  made  where  Vm ( t)  is  not  changing  rapidly.  The  total 
current  is  then 


I(t)  =  -V  (t)/(Ri  ||  R 2) 


(25) 


The  exact  equations,  including  the  time  constants,  are 


Vc(t)  =  V  +  V  (t)  +  —  /*  V(t’)  dt' 

S  0  m  R2C  2  o  m 


I  ( t)  =  -Vm(t)/R2  -  Vs(  t)/R  ! 


(26) 


+  ^‘t/RlCl  [v^ik  t 


This  leads  to  *10%  corrections  to  the  above  equations  for  our  typical  data. 


2.4 


BEAN  DIAGNOSTICS. 


Linear  processes,  such  as  attachment,  can  be  measured  from  the 
total  dose  and  the  integral  of  the  ion  chamber  current.  Recombination  is  a 
nonlinear  process  and  requires  for  its  determination  that  the  current  and 
dose  rate  be  known  at  a  particular  time.  Thus  it  is  very  important  in  this 
experiment  to  measure  accurately  both  the  dose  and  waveform. 

The  beam  diagnostics  used  during  this  experiment  were: 

1.  A  Rogowski  coil  that  meaured  the  derivative  of  the  beam  cur- 
rent,  1^,  at  the  entrance  to  the  ion  chamber.  When  numeri¬ 
cally  integrated  this  gives  the  dose  rate  as  a  function  of 
time,  y(t). 

2.  Faraday  cups  that  measured  Ib ( t) ,  or  y(t),  at  the  rear  of  the 
chamber. 

3.  Thin  Cu  foil  calorimeter  that  measured  total  dose  at  the 
rear  of  the  chamber. 


I 

I 


4.  TLD's  that  measured  total  dose  at  the  rear  of  the  chamber. 

The  last  three  of  the  diagnostics  were  located  in  the  backplate  of  the  cham¬ 
ber,  which  served  as  a  stopping  block  for  the  beam.  A  thin  foil  served  as  a 
window  so  that  the  diagnostics  could  be  in  the  room  air.  This  allowed  easy 
access  to  the  TLD's  without  disturbing  the  air  sample  and  removed  the 
requirement  for  a  vacuum  feedthrough  for  the  thermocouple  on  the  calorime¬ 
ter. 


2.4.1  Rogowski  Coil. 

The  Rogowski  coil  was  used  as  an  I  measuring  device  to  avoid  cir¬ 
cuit  oscillations  that  can  occur  when  the  usual  low  value  shunt  resistor  is 
used.  The  coil  consisted  of  a  simple  torus  of  teflon  insulated  wire  wound 
on  a  teflon  core.  The  area  of  a  turn  was  “0.4  cm2  and  there  were  about 
8  turns.  The  Rogowski  coil  was  mounted  in  a  housing  that  served  as  the 
final  length  of  the  collimator  in  front  of  the  ion  chamber.  The  diameter  of 
the  torus  was  1  3/8". 


It  can  be  shown  that  a  changing  current  passing  through  the  hole 
of  the  torus  produces  a  voltage  source  in  the  coil  that  is,  to  first  order, 
independent  of  position  of  the  current  and  is  given  by  KI  where 


K 


VA 

S 


(27) 


Here  n  is  the  number  of  turns,  A  is  the  area  of  a  turn,  and  S  is  the  circum¬ 
ference  of  the  torus.  The  equivalent  circuit  is  shown  in  Figure  6.  The 
inductance  of  the  Rogowski  coil  is 


m 


J 


r2  =  50  n 


Figure  6.  Rogowski  coil  equivalent  circuit. 


Solving  the  circuit  for  I,  we  find 


I(t)  = 


R,  +  R2  t  . 

_  /  V  (t’)dt'  +  — —  V  (t) 

KR,  o  m  KR,  m 


Using  the  values  given  above  we  find  that  K  =  3.7  nH,  L  =  30  nH  and 


I (t)  =  0.54  /  V  (t')dt’  +  0.16  V  (t) 


where  t  must  be  expressed  in  ns.  At  the  peak  of  the  pulse  the  integral  of 


V  is  at  least  10  times  larger  than  V  and  the  second  term  is  small  in  com- 
m  3  m 


pari  son  to  the  first.  The  second  term  represents  high  frequency  components 
that  are  lost  in  the  integration. 


The  Rogowski  coils  were  calibrated  by  running  a  square  pulse  of 
current  through  them  on  a  wire.  A  coaxial  ground,  short  length  of  wire,  and 
terminating  plate  with  50-n  resistor  were  needed  to  preserve  the  risetime  of 
the  pulse.  Figures  7  to  9  show  the  calibration  pulse,  the  Rogowski  coil 
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Figure  7.  Calibration  pulse  for  Rogowski  coil 


Numerical  integration  of  Rogowski  coil.  A  slight 
baseline  shift  in  the  Rogowski  coil  waveform  causes 
the  baseline  of  the  integrated  waveform  to  be  slanted 
this  was  taken  into  account  in  the  calibration 
procedure. 


waveform,  and  the  integral  of  the  Rogowski  coil.  Note  that  the  rise  and 
fall  times  of  the  reconstructed  pulse  are  very  close  to  those  of  the  input 
pulse  but  that  the  high  frequency  wiggles  on  the  top  of  the  pulse  are  not  as 
pronounced.  An  improvement  could  be  made  by  adding  in  part  of  the  un- inte¬ 
grated  waveform  as  in  the  second  term  in  Equation  (29),  but  it  was  not 
judged  to  be  necessary  for  the  purpose  of  this  experiment.  Very  accurate 
waveforms  of  the  HIFX  current  are  obtained  with  the  Faraday  cups.  Using 
just  the  integral,  the  coefficient  of  the  first  term  in  Equation  (29)  was 
found  to  be  0.36,  0.40,  0.41  A/V*ns  for  the  three,  nominally  identical 
Rogowski  coils  used  here. 

2.4.2  Faraday  Cup. 

In  a  corruption  of  the  original  meaning,  any  device  used  to  stop 
an  electron  beam  and  measure  the  current  is  nowadays  called  a  Faraday  cup. 
It  is  very  difficult  to  get  the  inductance  low  enough  to  avoid  oscillations 
in  a  Faraday  cup  used  to  stop  a  large  fraction  of  the  beam  and  shunt  it  to 
ground  through  a  low  impedance.  We  have  taken  the  approach  of  using  a  50- a 
cable  as  the  load  on  the  Faraday  cup  and  then  choosing  the  area  small  enough 
that  the  developed  voltage  is  below  the  breakdown  potential  of  the  cabling. 
The  50  n  impedance  insures  efficient  damping  of  the  leakage  flux  caused  by 
non-50  «  geometry  at  the  cup. 

The  Faraday  cups  used  here  were  simple  devices  fabricated  from 
0.141  Cujac  cable  and  small  A1  pellets  and  were  mounted  directly  in  a  hole 
in  the  back  plate  of  the  chamber  as  shown  in  Figure  10.  The  A1  pellet  was 
held  to  the  center  condcutor  of  the  Cujac  by  heating  the  pellet,  cooling  the 
wire,  and  allowing  both  to  come  to  room  temperature  after  insertion.  The 
pellet  was  completely  wrapped  in  teflon  tape  to  avoid  shorting  to  the  back- 
plate  or  foil  window.  The  pellet  is  more  than  an  electron  range  thick  and 
the  Faraday  cups  were  tested  to  withstand  500  V.  The  area  of  the  Faraday 
cups  was  typically  0.07  cm2;  a  100  A/cm2  current  at  the  rear  of  the  chamber 
would  produce  a  voltage  of  350  V.  During  these  experiments  voltages  up  to 
800  V  were  measured  without  any  evidence  for  breakdown. 


CONDUCTIVE 
„  EPOXY 


CUJAC 


Figure  10.  Faraday  cup  construction.  Twice  scale. 

A  calculation  showed  that  the  signal  due  to  electrons  stopping  in 
the  teflon  was  only  about  4X  of  that  due  to  electrons  collected  by  the  A1 
pellet.  This  is  offset  by  the  loss  of  charge  due  to  back  scattering  and 
secondaries  which  amounts  to  about  5%.  The  capacitance  of  the  Faraday  cup 
can  be  estimated  from  the  dimensions  of  Figure  10  to  be  about  4  pF. 

The  equivalent  circuit  of  the  Faraday  cup  is  shown  in  Figure  11. 


Figure  11.  Equivalent  circuit  of  Faraday  cup. 

The  current  source  of  the  incident  beam  causes  electrons  to  flow 
back  to  ground  through  the  stray  capacitance  and  terminating  resistor.  We 
find 
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where  t  =  0.2  ns.  The  Faraday  cup  measurements  were  normally  used  at  the 
peak  of  the  pulse  or  some  other  place  where  dV/dt  is  small.  The  second  term 
is  normally  2%  or  less  of  the  first  term  and  was  neglected  here. 


2.4.3  Calorimeter. 

The  5  mil  thick  Cu  foil  calorimeter  is  thin  compared  to  the  15  mil 
average  range  of  2  MeV  electrons  in  Cu.  At  2  MeV  the  stopping  power  is  a 
very  weak  function  of  the  energy.  A  chromel-al umel  thermocouple,  with  a 
sensitivity  of  40  uV/°C  at  room  temperature,  was  attached  to  the  foil  with  a 
small  drop  of  thermally  conducting  epoxy.  The  foil  had  a  diameter  of  6  mm, 
was  sandwiched  between  two  styrofoam  sheets,  and  mounted  in  a  hole  in  the 
backplate. 

The  mass  of  Cu  was  more  than  100  times  the  mass  of  the  epoxy  plus 
thermocouple  leads,  so  the  specific  heat  of  Cu  (0.092  cal/gm  °C)  controlled 
the  temperature  rise  of  the  calorimeter.  The  major  source  of  heat  loss  from 
the  calorimeter  is  down  the  thermocouple  leads  and  a  calculation  showed  a 
time  constant  of  6  s.  This  is  in  excellent  agreement  with  the  value  of  11  s 


found  from  the  semi -log  plot  of  a  calorimeter  record  in  Figure  12. 
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The  thermocouple  leads  ran  to  a  junction  box  about  2  ft  from  the 
beam  line  where  the  chromel  and  alumel  wires  were  connected  to  the  copper 
wires  of  a  twisted  shielded  pair  (TSP)  cable.  These  connections  were  made 
on  a  massive  Cu  block  to  insure  that  the  secondary  junctions  remained  at 
room  temperature.  Both  sides  of  the  circuit  were  also  capacitively  coupled 
to  ground  at  this  point  to  remove  any  fast  transient  introduced  by  the  beam. 
The  TSP  cable  connected  to  an  operational  amplifier  in  the  screen  room  and 
then  to  a  Visicorder. 
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The  dose  and  dose  rates  used  in  this  experiment  are  larger  than 
those  encountered  in  most  TLD  applications.  It  was  decided  to  use  a  TLD 
composed  of  5%  CaF2:Mn  powder  in  a  Teflon  matrix  because  there  was  more  data 
on  this  material  at  high  dose  and  high  dose  rate  (Ref.  1)  than  newer  mater¬ 
ials  that  may  have  superior  properties.  This  material  showed  no  dose  rate 
dependence  at  3  x  10 11  rad/s  from  Aurora  in  the  bremsstrahlung  mode  but 
starts  to  saturate  at  100  krads  total  dose. 

The  available  TLD's  were  6  mm  in  diameter  and  16  mm  thick.  In 
order  to  map  the  dose  at  the  rear  of  the  ion  chamber  the  TLD's  were  punched 
to  a  3  mm  diameter.  The  small  TLD's  were  recalibrated  by  HDL  personnel  with 
a  Co60  source  to  300  krads  and  were  found  to  scale  from  the  usual  TLD  cali¬ 
bration  curve  by  the  ratio  of  the  mass  of  the  small  TLD's  to  the  large 
TLD's.  The  TLD's  were  wrapped  in  A1  foil,  glued  to  an  A1  screw,  and  inserted 
into  the  backplate  from  the  rear.  The  TLD's  were  used  on  a  number  of  shots 
as  a  cross  check  of  the  other  beam  diagnostics,  but  the  shot- to-shot  mea¬ 
surement  of  the  dose  was  made  with  the  calorimeter,  Faraday  cups  and  Rogow- 
ski  coil. 

2.4.5  Layout. 

The  various  beam  diagnostics,  except  for  the  Rogowski  coil,  were 
mounted  in  holes  in  the  backplate  of  the  chamber  and  distributed  in  order  to 
check  the  uniformity  of  the  beam.  The  pattern  is  shown  in  Figure  13. 


T  =  TLD 

FC  =  FARADAY  CUP 
CAL  =  CALORIMETER 


Figure  13.  Pattern  of  beam  diagnostics  in  the  chamber  backplate. 

2.5  THE  EQUILIBRIUM  FIELO  MEASUREMENT. 

A  check  can  be  made  on  the  electron-ion  recombination  coefficient 
mesaured  with  the  ion  chamber  by  directly  generating  and  measuring  the  equi¬ 
librium  field.  The  basic  idea  of  the  way  that  this  can  be  done  is  shown  in 
Figure  14. 


Figure  14.  Simplified  diagram  of  measurement  of  the  equilibrium  field. 


This  device,  which  we  refer  to  as  an  equilibrium  field  diode,  consists  of 
two  electrodes:  a  thin  grounded  electrode  and  a  thick  high  voltage  elec¬ 
trode.  The  electron  beam  passes  through  the  front  electrode  but  is  stopped 
by  the  rear  electrode.  The  electron  beam  serves  not  only  to  ionize  the  air 
but  also  to  produce  a  field  by  charging  the  thick  electrode.  As  the  thick 
electrode  charges  and  a  field  builds  up  in  the  air,  a  conduction  current 
carries  charge  through  the  air  back  to  the  grounded  thin  electrode. 

The  potential  of  the  thick  electrode  approaches  the  equilibrium 
voltage  according  to 

VS  «  RAIb  (1  -  e"t/RAC)  (31) 

where  is  the  resistance  of  the  air  and  C  is  the  capacitance  between  thick 
and  thin  electrodes.  If  the  dimensions  are  similar  to  those  of  the  ion 
chamber,  the  air  resistance  will  be  »Z  a  at  the  highest  dose  rate  and  C  « 
10  pF.  This  shows  that  the  current  drawn  through  the  measuring  voltage 
divider  will  be  negligible  compared  to  the  current  through  the  air  and  that 
the  steady  state  voltage,  Vs>  is  attained  with  a  very  short  time  constant, 
«0.2  ns. 


As  with  the  ion  chamber  a  boundary  layer  will  form  at  the  thick 
electrode  and  its  potential  will  be  larger  than  that  given  by  Equation  (31). 
This  can  again  be  avoided  by  measuring  the  voltage  to  the  thin  electrode 
from  a  grid  inserted  between  the  two  electrodes. 

2.5.1  Matching  «c. 

The  equilibrium  field  diode,  as  described  so  far,  would  not  mea¬ 
sure  the  equilibrium  field  generated  by  a  nuclear  burst.  This  is  because 
the  Vilue  of  Kc  in  Equation  (4)  is  larger  for  a  parallel  beam  of  electrons 
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than  for  the  divergent  flux  of  Compton  electrons  produced  by  gamma  rays. 

Since  the  air  would  be  receiving  a  smaller  dose  rate  for  the  same  current 

density  the  fields  in  the  equilibrium  field  diode  would  be  larger.  For 

y  rays,  K  =  2.5  x  10" 12  (A/cm2)/(rad/s) ,  whereas  for  HIFX,  K  «5.9  x 
c  c 

10-12  ( A/cm2) /( rad/s).  The  equilibrium  fields  from  gammas  and  from  an  elec¬ 
tron  beam  differ  by  even  more  than  the  ratio  of  the  K  ‘s  because  of  the 

c 

electric  field  dependence  of  the  electron  mobility.  The  equilibrium  field 
diode  was  compensated  in  this  experiment  for  this  effect  by  making  the  thick 
electrode  transmitting  over  part  of  its  area  and  adding  a  second  ground 
electrode  as  in  Figure  15. 

A  sector  hole  was  cut  in  the  center  electrode,  but  the  hole  is 
covered  on  both  sides  by  1  mil  A1  foils.  Only  part  of  the  beam  now  acts  to 
charge  the  thick  electrode,  and  an  additional  air  resistance  is  introduced 
in  the  rear  of  the  diode.  The  size  of  the  hole  is  chosen  to  make  the  field 
the  same  as  that  from  gamma  rays.  The  steady  state  electric  field  is  given 
by 


Es  =  Wd 


(32) 


where  d  is  the  spacing  between  electrodes.  Allowing  for  an  area  of  beam 
stopping  ( A t )  different  from  the  area  of  the  resistor  (A2)  we  have 


E  =  Jb  A^ 


(33) 


We  call  the  fraction  of  the  area  of  the  thick  electrode  that  is  transparent 
to  the  beam  n.  Then  A  =  (1  -  n)  A  and  A  =  (1  +  n)A  ,  where  A  is  the 
area  of  the  beam.  Equation  (33)  becomes 
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Since  for  HIFX  electrons  is  2.4  times  too  large,  we  chose  n  =  0.42  (angle 
of  sector  hole  =  150°)  to  make  the  diode  produce  the  same  fields  that  would 
be  produced  by  gamma  rays. 

The  whole  equilibrium  field  diode  assembly  was  covered  by  a  cylin¬ 
der  rolled  from  a  sheet  of  Cu  foil  with  holes  for  the  center  electrode  and 
grid  leads.  This  cylinder  served  as  the  return  current  path  for  the  beam, 
keeping  any  back  emf  small  compared  to  the  voltage  developed  in  the  diode. 
The  diode  was  mounted  in  place  of  the  ion  chamber  in  the  vacuum  chamber 
(See  Figure  2).  With  the  center  electrode  partially  transmitting,  the  beam 
diagnostics  in  the  backplate  could  be  used. 


2.5.2 


Elevated  Gas  Temperature. 


The  recombination  process  with  water  molecules  acting  as  a  third 
body  depends  upon  the  loss  of  energy  by  the  electron  to  the  rotational 
levels  of  the  H20  molecule.  We  have  seen  that  there  is  theoretical  and 
experimental  evidence  to  suspect  that  the  value  of  depends  strongly  on 
the  temperature  of  the  electrons.  A  somewhat  heuristic  extention  of  the 
Thomson  theory  of  three-body  recombination  (Ref.  6)  led  to  the  expression 


g3  =  5.7  x  10-9  Qhj  T-3  (cm6/s) 


(35) 


where  v  /N  is  the  energy  exchange  rate  coefficient  and  can  be  found  from 
electron  swarm  data.  In  this  derivation  the  electron  temperature  is  taken 
as  the  gas  temperature,  T. 
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Since  the  spacing  of  the  rotation  energy  levels  in  asymmetric  top 
molecules  like  H2O  are  on  the  order  of  10" 3  to  10-4  eV  (Ref.  3),  the  distri¬ 
bution  of  the  H20  molecules  among  the  rotational  levels  can  be  changed  by 
modest  changes  in  gas  temperature.  Thus  there  is  a  possibility  that  besides 
the  dependence  of  63  on  electron  temperature  there  is  also  a  dependence  on 
neutral  gas  temperature  through  the  v^/N  factor  in  Equation  (35).  The  equi¬ 
librium  field  diode  was  designed  to  operate  at  room  temperature  and  at  a 
temperature  of  200°C.  At  the  fields  generated  in  the  diode,  the  electron 
temperature  is  determined  by  the  field,  and  a  large  change  in  the  steady 
state  voltage  can  be  attributed  to  the  dependence  of  B3  on  gas  temperature. 
Due  account  must  be  taken  for  other  temperature  dependent  quantities  such  as 
the  electron  mobility. 

A  thermistor  was  epoxied  into  a  hole  in  the  Rogowski  coil  housing 
and  monitored  in  the  screen  room  during  high  temperature  operation  of  the 
equilibrium  field  diode. 


SECTION  3 
DATA  ANALYSIS 


3.1  DOSIMETRY. 

Two  configurations  were  used  to  transport  the  electron  beam  into 
the  chamber.  In  the  first  arrangement  the  vacuum  chamber  was  mounted  on  the 
end  of  the  drift  tube  and  only  a  2  mil  A1  window  separated  the  drift  tube 
partial  pressure  from  the  atmospheric  pressure  of  the  chamber.  The  dose 
calculated  from  an  integration  of  the  Faraday  cup,  the  calorimeter,  and  the 
TLD’s  agreed  to  within  10%.  The  dose  obtained  from  a  double  integration  of 
the  Rogowski  coil  was  a  factor  of  40%  larger,  which  is  consistent  with  a 
calculation  of  the  scattering  of  the  electron  beam  in  the  foils  and  grids  in 
the  chamber. 

In  the  second  configuration,  there  was  a  3-cm  air  gap  between  the 
end  of  the  beam  tube  and  the  vacuum  chamber.  The  calorimeter  now  measured 
40%  larger  than  the  central  Faraday  cup,  and  the  Rogowski  coil  was  found  to 
give  a  dose  1.82  ±  0.15  times  that  of  the  calorimeter  on  27  shots  for  which 
the  Rogowski  coil  signal  was  manually  digitized  and  doubly  integrated.  This 
second  configuration  is  the  one  that  was  used  for  the  high  dose-rate  data 
from  which  the  electron-ion  recombination  coefficient  was  found. 

In  analyzing  the  data  the  chamber  was  divided  into  a  front  half 
with  the  dose  rate  given  by  the  Rogowski  coil  and  a  back  half  with  the  dose 
rate  scaled  to  the  calorimeter.  The  dose-times-area  product  was  kept  con¬ 
stant  in  keeping  with  the  assumption  that  beam  spreading  was  responsible  for 
the  fall-off  in  the  dose. 
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3.2  GRID  SIGNALS. 

The  voltage  signal  measured  by  the  grid  has  many  potential  sources: 


» 


1. 


I 

I 
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2. 


i 


i 


3. 


4. 


The  applied  voltage. 

Drop  in  applied  voltage  to  the  air  as  the  chamber  draws  cur¬ 
rent  through  the  chamber  load  resistance. 

Magnetic  coupling  to  the  beam. 

Charge  in  transit  through  the  chamber. 


i 


The  first  three  of  these  have  already  been  treated. 


The  beam  that  traverses  the  chamber  is  neutralized  by  plasma  elec¬ 
trons  that  are  driven  to  the  walls,  however,  there  is  a  time  constant  of 
t  =  eQ/a  needed  to  relax  charge  injected  into  the  chamber  by  a  change  in  the 
beam  current.  For  the  initial  rise  of  the  current  there  is  an  additional 
time  required  to  build  up  the  conductivity  in  the  air.  Once  there  is  suffi¬ 
cient  conductivity,  the  rate  equation  for  the  net  charge  in  the  chamber  is 

dp  =  Il£  dt  -  7«oE  dt  (36) 

c  dt 

i 

i 


. 

1 


where  J/c  is  the  charge  density  in  the  beam.  Rearranging, 


If  t  is  much  less  than  the  time  interval  over  which  J  changes,  the  quasi¬ 
steady  state  solution  for  p  is 


J 

m 

m 

I 

i 


1 


Since  the  grid  will  see  a  potential  proportional  to  p,  it  will  have  a  compo¬ 
nent  that  goes  as  the  derivative  of  the  beam  current. 

In  Figure  16a  is  shown  a  grid  voltage  signal  produced  on  a  shot 
with  Vq=  0.  The  peaks  and  valleys  can  be  associated  with  similar  features 
in  the  derivative  of  the  beam  current  given  by  Rogowski  coil  signal  of  Fig¬ 
ure  16b.  In  the  example  following  Equation  (14)  we  estimated  t  to  be  5  ps 
and  found  that  a  field  on  the  order  of  100  V/cm  would  be  produced.  This  is 
of  the  same  size  as  the  field  measured  in  Figure  16a. 

The  largest  uncertainty  in  this  experiment  is  the  voltage  on  the 
air  as  measured  by  the  grids.  The  quality  of  the  grid  signal  increased  as 
the  applied  voltage  increased.  This  is  because  the  component  of  grid  volt¬ 
age  due  to  charge  injected  into  the  chamber  decreases  in  relative  importance 
and  in  absolute  value  as  the  air  conductivity  is  increased  at  higher  field. 

Figure  17  shows  that  at  an  applied  voltage  of  «4  kV,  the  grid  voltage  could 

be  determined  quite  accurately. 

3.3  STOPPING  CURRENT. 

Charge  from  the  beam  that  stops  in  the  air,  grids,  and  foils  of 
the  chamber  will  produce  signals  even  when  the  applied  voltage  is  zero  as 
shown  by  Figure  18.  These  stopping  currents  carry  charge  from  ground  to 

grid  and  center  electrode  and  can  be  included  in  the  circuit  as  in  the  dia- 

S  S 

gram  of  Figure  19.  Ic  E  is  the  stopping  current  to  the  center  electrode,  IG 

is  the  stopping  current  to  the  grid,  is  the  resistance  of  the  air,  and 

Rg  L  is  the  resistance  of  the  boundary  layer.  The  node  equations  lead  to 


hyfnft 


GRID  VOLTAGE  (85  V/div) 


SHOT  345 


V  =  0 
o 


y(peak)  =  6.7x10  rad/s 


ION  CHAMBER  CURRENT 


GRID  VOLTAGE 


TIME  (10  ns/div) 


Figure  18.  Signals  produced  by  beam  stopping  currents. 
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GRID  VOLTAGE  (85  V/div) 


Assume 


I 


S 

C.E. 


=  K 


S- 

Y 


Then  for  a  VQ 


=  0  shot  we  can  find  K 


S 


from 


K  *  '-VRL  -  VV^ 


Once  we  know  K  the  corrected  current  through  the  air  is 


*A  "  <W  + 


Figure  19.  Ion  chamber  circuit  with  stopping  currents  included 


One  problem  in  finding  K  from  Equation  (40)  is  that  it  requires 
the  air  resistance,  which  is  not  know  a  priori.  After  an  initial  guess  at 
R^,  the  data  was  analyzed  and  a  value  for  was  found.  The  data  was  then 
analyzed  a  second  time  using  the  value  of  in  Equation  (40)  from  the  first 
value  found  for  The  grid  voltage  trace  shown  in  Figure  18  is  the  best 
voltage  record  obtained  on  a  V0  =  0  shot  (because  of  the  smooth  beam  cur¬ 
rent  waveform).  Using  this  voltage,  the  recorded  value  of  y,  and  the  field 
dependent  values  of  u,  a,  and  ge,  we  find  an  air  resistance  of  “3.5  n  and 

c 

from  Equation  (40)  a  value  of  4.6  x  10" 13  A/( rad/s)  for  K  .  The  ion  chamber 
current  was  then  corrected  according  to 

!a  =  VRL  +  4‘6  x  10-13  ^  (42) 

where  y  is  measured  by  the  Rogowski  coil. 

If  one  uses  the  rule  of  thumb  that  the  fraction  of  beam  electrons 

stopped  in  a  material  is  equal  to  the  material's  thickness  expressed  as  a 

c 

fraction  of  the  electron  range,  the  value  of  K  estimated  for  the  grid,  cen¬ 
ter  electrode  foil,  and  air  in  the  chamber  is  x  10- 13  A/(rad/s),  in  good 
agreement  with  the  value  found  above. 

3.4  THE  ELECTRON-ION  RECOMBINATION  COEFFICIENT. 

We  assume  that  the  electron  density  in  the  ion  chamber  is  in  a 

quasi-steady  state,  that  is,  that  the  characteristic  time  for  a  change  in 

dose  rate  is  long  compared  to  the  relaxation  time  of  the  electron  density 

(t  =4 /a  +  1 / 6e n ) .  If  the  dose  is  uniform  throughout  a  volume, 

n  =  I/euEA  (43) 


>yiviVi»,3 

'  >■) 


2(K y  -  an) 


(44) 


0e  = 


n2  /l  +  /l  +  4  a  \ 

l  "V 


In  order  to  take  into  account  two  regions  with  different  dose 
rates  we  used  the  following  procedure: 


1.  Guess  nt. 

2.  Calculate  r\2  =  1  /A 2( I /e yE  -  n1A1/A2). 

3.  Generate  B1  and  B2  from  Equation  (44). 

e  e 

4.  Iterate  until  B1  =  32. 

e  e 

Previously  measured,  field-dependent  values  of  w  and  a  were  used  (Ref.  7). 


The  graph  of  Figure  20  shows  the  resulting  electron- ion  recombina¬ 
tion  coefficient  as  a  function  of  electric  field  over  the  range  that  could 
be  attained  in  the  ion  chamber.  The  moist  air  data  appears  to  be  approach¬ 
ing  the  expected  three-body  value  { «1 . 7  x  10-5  cm3/s)  at  low  field.  At 
fields  of  interest  to  EMP  (4000  V/cm),  however,  has  fallen  to  4  x 
10-6  cm3/s,  a  value  consistent  with  clustered  two-body  recombination,  only. 
The  dry  air  data  are  about  a  factor  of  ten  below  the  moist  air  at  300  V/cm 
where  the  data  overlap.  The  dry  air  data  spans  a  range  of  electric  field 
lower  than  that  of  the  moist  air.  This  is  because  the  decreased  values  of 
Be  correspond  to  larger  chamber  currents  and  thus  more  voltage  dropped  e 
across  the  chamber  load  resistor. 


The  error  bars  in  Figure  20  were  based  on  fixed  uncertainties  of 
10%  for  u,  a,  and  y»  5%  for  A,  and  100%  for  B^ .  However,  the  error  bars  are 
really  determined  by  the  uncertainties  in  the  chamber  current  and  voltage. 
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The  uncertainty  in  the  current  was  taken  to  be  10%  of  the  uncorrected  cur¬ 
rent  plus  30%  of  the  correction.  The  uncertainty  in  the  voltage  varied  from 
shot  to  shot  depending  on  the  shape  of  the  individual  dose-rate  waveform  and 
applied  voltage;  values  varied  from  10  to  100%. 

Some  data  were  taken  in  addition  to  that  shown  in  Figure  20.  The 
ion  chamber  with  the  2  cm  spacing  was  used  extensively  in  trying  to  get  data 
at  lower  fields.  Unfortunately,  the  divergence  of  the  electron  beam  was 
such  that  some  beam  electrons  stopped  in  the  teflon  spacers,  making  very 
large  zero  voltage  signals.  Some  data  were  taken  at  lower  dose  rates, 
«10 1 2  rad/s.  These  data  verified  the  expected  value  of  a  to  within  about 
10%,  but  it  is  not  possible  to  obtain  accurate  recombination  information, 
because  recombination  is  responsible  only  for  about  10%  of  the  free  electron 
removal . 

3.5  THE  EQUILIBRIUM  FIELD. 

The  equilibrium  field  diode  was  also  exposed  in  the  two  configura¬ 
tions:  when  mounted  on  the  end  of  the  60-cm  drift  tube  the  entrance  dose 
rate  was  «2  x  10 12  rad/s,  when  mounted  separate  from  the  10  cm  drift  tube 
the  dose  rate  was  «5  x  10 13  rad/s.  The  grid  voltage  was  again  very  uncer¬ 
tain  because  of  noise  due  to  charge  injected  into  the  chamber.  One  of  the 
best  grid  voltage  records  is  shown  in  Figure  21a;  a  much  more  uncertain 
voltage  signal  is  shown  in  the  second  signal  in  Figure  21b. 

Table  1  lists  the  values  of  the  field  for  the  high  and  low  expo¬ 
sure,  for  moist  and  dry  air  and  for  room  temperature  and  200°C.  We  used  the 

values  of  e  measured  with  the  ion  chamber  with  two  different  dose-rate 
e 

regions  (front  and  back)  and  predicted  the  equilibrium  field.  An  iterative 
procedure,  like  that  for  the  ion  chamber  analysis,  was  needed  because  of  the 
field  dependent  quantitites  (a,  u,  f3g)  involved  in  calculating  Es»  For  the 


lower  dose  rate  the  field  is  determined  mostly  by  attachment.  Our  predicted 
value  is  950  V/cm,  about  20 %  higher  than  the  measured  values  in  Table  1.  At 
the  high  dose  rate  the  field  is  determined  mostly  by  recombination.  Our 
predicted  value  is  2300  V/cm,  about  50%  higher  than  the  measured  value.  If 
(because  of  three-body  recombination)  the  value  of  was  2  x  10* 5  cm3/s  the 
predicted  equilibrium  field  at  the  high  dose  would  be  greater  than  7  kV/cm. 

Table  1.  The  Measured  equilibria  fields. 


Shot 

• 

Y 

(rad/s) 

E 

V/cm 

GAS 

274 

»1  x  io12 

840 

2.4%  H2O 

275 

2.2  x  1012 

770 

II 

276 

2.1  x  1012 

740 

II 

277 

2.0  x  1012 

770 

II 

283 

4.5  x  1013 

1700 

II 

286 

5.3  x  1013 

“1800 

II 

287 

4.4  x  1013 

1300 

II 

288 

4.4  x  1013 

“2000 

2.1%  H20,  210°C 

289 

4.9  x  1013 

II 

293 

4.6  x  1013 

Dry  air,  200°C 

The  high  temperature  measurements  were  very  uncertain  and  it  is 
not  possible  to  tell  if  there  was  a  change  in  the  equilibrium  field.  Since 
the  temperature  effect  of  interest  involves  three-body  recombination  and  we 
have  convincing  evidence  that  there  is  little  of  the  three-body  process 
operating  at  these  fields,  we  would  not  expect  to  see  large  changes  in  the 
equilibrium  field. 


SECTION  4 
CONCLUSIONS 


This  experiment  shows  that,  with  careful  design,  large  ionization 
can  be  produced  and  measured  in  an  ion  chamber  using  a  high  current  electron 
beam.  Magnetic  coupling  can  be  kept  to  a  minimum  by  keeping  the  beam  diame¬ 
ter  small  and  the  current  return  path  close.  A  field  proportional  to  the 
derivative  of  the  beam  current  due  to  the  finite  relaxation  time  of  the 
charge  injected  into  the  chamber  is  the  largest  perturbation  caused  by  the 
electron  beam.  This  makes  it  important  to  try  to  keep  the  current  waveform 
smooth  without  an  excessively  fast  risetime. 

It  was  also  shown  that  very  large  current  densities  can  be  mea¬ 
sured  with  an  ion  chamber.  The  boundary  layer  did  not  consume  the  whole 
applied  voltage.  A  significant  amount  of  the  voltage  appeared  on  the  boun¬ 
dary  layer,  but  sufficient  fields  remained  in  the  bulk  of  the  gas.  Boundary 
layer  instabilities  did  not  result  in  breakdown  occurring  between  the 
cathode  and  anode  of  the  chamber. 

Both  the  ion  chamber  measurement  of  the  electron-ion  recombination 
coefficient  and  the  measurement  of  the  equilibrium  field  show  that  three- 
body  recombination  is  not  a  significant  electron  removal  process  at  EMP 
fields. 
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NAVAL  AIR  SYSTEMS  COMMAND 
ATTN:  AIR350F 
ATTN:  AIR  5161 

NAVAL  FACILITIES  ENGINEERING  COMMAND 
ATTN:  04E 
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NAVAL  INTELLIGENCE  SUPPORT  CTR 
ATTN:  NISC-45 


NAVAL  OCEAN  SYSTEMS  CENTER 

ATTN:  COOE  4471  TECH  LIB 


NAVAL  ORDNANCE  STATION 

ATTN:  STANDARDIZATION  DIVISION 


NAVAL  POSTGRADUATE  SCHOOL 

ATTN:  COOE  1424  LIBRARY 


NAVAL  RESEARCH  LABORATORY 

ATTN:  CODE  2627  TECH  LIB 


NAVAL  SURFACE  WEAPONS  CENTER 
ATTN:  CODE  F30 
ATTN:  COOEF32 
ATTN.  CODE  H20 
ATTN.  CODE  R40 
ATTN:  CODE  R43 
ATTN:  CODE  425 


NAVAL  SURFACE  WEAPONS  CENTER 
ATTN:  CODE  F  56 


NAVAL  TELECOMMUNICATIONS  COMMAND 
ATTN:  DEPUTY  DIRECTOR  SYSTEMS 


NAVAL  WEAPONS  CENTER 

ATTN:  CODE  343  FKA6A2  TECH  SVCS 


OFC  OF  THE  DEPUTY  CHIEF  OF  NAVAL  OPS 

ATTN.  NOP  098  OFC  RES- DEV  TEST  &  EVAL 
ATTN:  NOP  506 
ATTN:  NOP  551 

ATTN:  NOP  654  STRAT  EVAL  &  ANAL  BR 
ATTN:  NOP  94 
ATTN.  NOP  981 
ATTN:  N0P981N1 


OFFICE  OF  NAVAL  RESEARCH 
ATTN.  CODE  427 


SPACE  &  NAVAL  WARFARE  SYSTEMS  CMD 
ATTN:  PME  117-21 
ATTN:  TECHNICAL  LIBRARY 


STRATEGIC  SYSTEMS  PROGRAM S(PM- 1) 
ATTN:  NSP  43  TECH  LIB 


THEATER  NUCLEAR  WARFARE  PROGRAM  OFC 
ATTN:  PMS423 


U  S  NAVAL  FORCES,  EUROPE 
ATTN:  N54 


DEPARTMENT  OF  THE  AIR  FORCE 


AERONAUTICAL  SYSTEMS  OIVISION,  AFSC 
ATTN:  ASD/ENESS  P  MARTH 
ATTN:  ASD/ENSSA 
ATTN:  ASO/YYEF 


AEROSPACE  DEFENSE  COMMAND 
ATTN:  ADC0M/J2A 
ATTN:  J6T 


AIR  FORCE/INE 

ATTN:  INA 


AIR  FORCE  AERONAUTICAL  SYS  DIV 
ATTN:  AFWAL/FIEA 
ATTN:  ASO/ENACE 


AIR  FORCE  COMMUNICATIONS  COMMAND 
ATTN:  SIM/SA 


AIR  FORCE  GEOPHYSICS  LABORATORY 
ATTN:  SULL 


AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 
ATTN:  LIBRARY 


AIR  FORCE  WEAPONS  LABORATORY.  AFSC 
ATTN:  NAT 
ATTN:  NT 
ATTN:  NTAEE 
ATTN:  NTYC 
ATTN:  NTYE 
ATTN:  NTYEE 
ATTN:  NTYEP 
ATTN:  SUL 


AIR  LOGISTICS  COMMAND 
ATTN:  OO-ALC/MM 
ATTN:  OO-ALC/MMETH 
ATTN:  OO-ALC/MMETH 


AIR  UNIVERSITY  LIBRARY 
ATTN:  AUL-LSE 


ASSISTANT  CHIEF  OF  STAFF 

ATTN:  AF/SAMI  TECH  INFO  DIV 


BALLISTIC  MISSILE  OFFICE/DAA 
ATTN:  ENSN 


DEPUTY  CHIEF  OF  STAFF 
ATTN:  LEEEU 


DEPUTY  CHIEF  OF  STAFF 
ATTN:  AF/RDQI 


DEPUTY  CHIEF  OF  STAFF 

ATTN:  AFRDS  SPACE  SYS  *  C3  DIR 


DEPUTY  CHIEF  OF  STAFF 

ATTN:  AFXOXFM  PLNS,  FRC  DEV  MUN  PLNS 


ELECTRONIC  SYSTEMS  DIVISION 
ATTN:  SCS-1E 


FOREIGN  TECHNOLOGY  OIVISION.  AFSC 
ATTN:  NIIS  LIBRARY 


UNITED  STATES  AIR  FORCE  HEADQUARTERS 
ATTN  SAS 


DEPARTMENT  OF  THE  AIR  FORCE  (CONTINUED) 

NORAD 

ATTN:  NORAD/J5YX 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
ATTN:  TSLD 

SPACE  COMMAND 
ATTN  DEE 

SPACE  DIVISION 

ATTN:  IND 

SPACE  DIVISION 

ATTN:  YKM 

STRATEGIC  AIR  COMMAND 
ATTN:  DEMUE 

STRATEGIC  AIR  COMMAND 
ATTN:  DEPR 

STRATEGIC  AIR  COMMAND 
ATTN:  DOCSD 

STRATEGIC  AIR  COMMAND 
ATTN:  INAD 

STRATEGIC  AIR  COMMAND 
ATTN:  INAO 

STRATEGIC  AIR  COMMAND 
ATTN:  SAC/ LG W 

STRATEGIC  AIR  COMMAND 
ATTN:  NRI/STINFO 

STRATEGIC  AIR  COMMAND 
ATTN:  SAC/SIP 

STRATEGIC  AIR  COMMAND 
ATTN:  XPFC 

STRATEGIC  AIR  COMMAND 
ATTN:  XPFR 

STRATEGIC  AIR  COMMAND 
ATTN:  XPFS 

STRATEGIC  AIR  COMMAND 
ATTN  XPQ 

TACTICAL  AIR  COMMAND 
ATTN  TAC/XPJ 

U  S  AIR  FORCES  IN  EUROPE 
ATTN  USAFE/XPXX 

U  S  RESEARCH  A  DEVELOPMENT  COORD 
ATTN  USRADCO 


DEPARTMENT  OF  ENERGY 

EMERGENCY  ELECTRIC  POWER  ADM 
ATTN:  LIBRARY 

UNIVERSITY  OF  CALIFORNIA 

LAWRENCE  LIVERMORE  NATIONAL  LAB 

ATTN:  L  658  TECH  INFO  DEPT  LIB 

LOS  ALAMOS  NATIONAL  LABORATORY 

ATTN:  MS  P364  REPORTS  LIBRARY 

SANDIA  NATIONAL  LABORATORIES 
ATTN:  TECH  LIB  3141 

OTHER  GOVERNMENT 

CENTRAL  INTELLIGENCE  AGENCY 
ATTN:  OSR/SE/C 
ATTN:  OSR/SE/F 
ATTN:  OSWR/NED 
ATTN:  OSWR/STD/MTB 

FEDERAL  EMERGENCYS  MANAGEMENT  AGENCY 
ATTN:  SL-EM 

NATIONAL  BUREAU  OF  STANDARDS 
ATTN:  723.03 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORP 

ATTN:  LIBRARY 

AGBABIAN  ASSOCIATES.  INC 
ATTN:  LIBRARY 

ALLIED  CORP 

ATTN:  DOCUMENT  CONTROL 

ALLIED  CORP 

ATTN:  DEPT  6401 

AT&T  TECHNOLOGIES.  INC 
ATTN  W  EDWARDS 

AVCO  SYSTEMS  DIVISION 

ATTN  LIBRARY  A830 

BDM  CORP 

ATTN  CORPORATE  LIB 

BDM  CORP 

ATTN  LIBRARY 

BOEING  CO 

ATTN  D  EGELKROUT 

ATTN  DKEMLE 

ATTN  H  WICKLEIN 

ATTN  KENT  TECHNICAL  LIBRARY 

ATTN  ORG  2  3744 

ATTN  R  SCHEPPE 
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INSTITUTE  FOR  DEFENSE  ANALYSES 
ATTN:  CLASSIFIED  LIBRARY 
ATTN:  TECH  INFO  SERVICES 


BOEING  MILITARY  AIRPLANE  CO 
ATTN:  C  SUTTER 

BOOZ-ALLEN  &  HAMILTON,  INC 

ATTN:  TECHNICAL  LIBRARY 

BOOZ-ALLEN  &  HAMILTON,  INC 
ATTN:  L  ALBRIGHT 

BOOZ  ALLEN  &  HAMILTON,  INC 
ATTN:  DDURGIN 

CALSPAN  CORP 

ATTN:  LIBRARY 

COMPUTER  SCIENCES  CORP 
ATTN:  ASCHIFF 

DIKEWOOD  CORP 

ATTN:  TECHNICAL  LIBRARY 

DIKEWOOD  CORP 

ATTN:  KLEE 


IRT  CORP 

ATTN:  B  WILLIAMS 
ATTN:  RW  STEWART 

ITT  TELECOMMUNICATIONS  CORP 
ATTN:  R  SCHWALLIE 

JAYCOR 

ATTN:  EWENAAS 

JAYCOR 

ATTN:  LIBRARY 

JOHNS  HOPKINS  UNIVERSITY 
ATTN:  P  PARTRIDGE 

KAMAN  SCIENCES  CORP 
ATTN:  LIBRARY 

KAMAN  SCIENCES  CORP 
ATTN:  E  CONRAD 


E  SYSTEMS,  INC 

ATTN:  E  WILKES 


KAMAN  TEMPO 

ATTN:  DASIAC 
ATTN:  R  RUTHERFORD 


E  SYSTEMS.  INC 

ATTN:  J  MOORE 


KAMAN  TEMPO 

ATTN:  DASIAC 


EG&G  WASH  ANALYTICAL  SVCS  CTR,  INC 
ATTN:  A  BONHAM 
ATTN:  C  GILES 


LITTON  SYSTEMS.  INC 

ATTN:  EEUSTIS 


ELECTRO  MAGNETIC  APPLICATIONS,  INC 
ATTN:  D  MEREWETHER 


LITTON  SYSTEMS,  INC 
ATTN:  J  MOYER 


FORD  AEROSPACE  &  COMMUNICATIONS  CORP 
ATTN:  H  LINDER 


LITTON  SYSTEMS,  INC 

ATTN:  J  SKAGGS 


GENERAL  ELECTRIC  CO 

ATTN:  C  HEWISON 


LOCKHEED  MISSILES  &  SPACE  CO,  INC 

ATTN:  TECH  INFO  CTR  D/COLL,  D/90  11 ,  B/ 106 


GRUMMAN  AEROSPACE  CORP 
ATTN:  L0135 


LTV  AEROSPACE  &  DEFENSE  COMPANY 
ATTN:  2  58010  LIBRARY 


HARRIS  CORP 

ATTN:  V  PRES  A  MGR  PRGMS  DIV 


LUTECH.  INC 

ATTN:  FTESCHE 


HERCULES,  INC 

ATTN:  W  WOODRUFF 


MCDONNELL  DOUGLAS  CORP 

ATTN:  TECHNICAL  LIBRARY  SERVICES 


HONEYWELL,  INC 

ATTN:  S&RC  LIBRARY 


MCDONNELL  DOUGLAS  CORP 

ATTN:  TECHNICAL  LIBRARY 


HONEYWELL,  INC 

ATTN:  LIBRARY 


METATECH  CORP 

ATTN:  R  SCHAEFER 


HUGHES  AIRCRAFT  CO 
ATTN:  CTD  6 


METATECH  CORPORATION 
ATTN:  WRADASKY 


IIT  RESEARCH  INSTITUTE 
ATTN:  IMINDEL 


MISSION  RESEARCH  CORP 
ATTN:  EMP  GROUP 
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MISSION  RESEARCH  CORP 
ATTN:  JLUBELL 

MISSION  RESEARCH  CORP.  SAN  DIEGO 
2  CYS  ATTN:  J  CHERVENAK 
2  CYS  ATTN:  V  VAN  LINT 

MITRE  CORP 

ATTN:  M  FITZGERALD 

PACIFIC  SIERRA  RESEARCH  CORP 

ATTN:  H  BRODE,  CHAIRMAN  SAGE 

PHYSICS  INTERNATIONAL  CO 

ATTN:  DOCUMENT  CONTROL 


RAD  ASSOCIATES 

ATTN  DOCUMENT  CONTROL 
ATTN:  PHAAS 
ATTN:  WKARZAS 

RAD  ASSOCIATES 

ATTN:  LIBRARY 

RAYTHEON  CO 

ATTN:  H  FLESCHER 

RCA  CORP 

ATTN:  GBRUCKER 

RESEARCH  TRIANGLE  INSTITUTE 
ATTN:  M  SIMONS 

ROCKWELL  INTERNATIONAL  CORP 

ATTN:  G  MORGAN  D256/MC  BA36 
ATTN:  J  ERB  D257/MC  BB17 

ROCKWELL  INTERNATIONAL  CORP 
ATTN:  B1DIVTICBAOB 


SCUBED 

ATTN:  A  WILSON 

SCIENCE  A  ENGRG  ASSOCIATES,  INC 
ATTN:  V  JONES 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  W  CHADSEY 

SCIENCE  APPLICATIONS.  INC 
ATTN:  EODONNELL 
ATTN:  TECH  INFO  CTR 

SPERRY  CORP 

ATTN.  R  LAZARCHIK 

SPERRY  CORP 

ATTN:  JINDA 

SPERRY  CORP 

ATTN:  TECHNICAL  LIBRARY 

SRI  INTERNATIONAL 

ATTN:  A  PADGETT 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  F  LEOPARD 

TEXAS  INSTRUMENTS.  INC 

ATTN:  TECHNICAL  LIBRARY 

TRW  ELECTRONICS  &  DEFENSE  SECTOR 
ATTN:  J  BROSSIER 
ATTN:  JPENAR 

TRW  ELECTRONICS  &  DEFENSE  SECTOR 
ATTN:  R  HENDRICKSON 

TRW  ELECTRONICS  &  DEFENSE  SECTOR 
ATTN:  LIBRARIAN 


